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ABSTRACT 

In a recent paper we have calculated the power density spectrum of Gamma-Ray Bursts arising from multiple 
shocks in a relativistic wind. The wind optical thickness is one of the factors to which the power spectrum is most 
sensitive, therefore we have further developed our model by taking into account the photon down-scattering on 
the cold electrons in the wind. For an almost optically thick wind we identify a combination of ejection features 
and wind parameters that yield bursts with an average power spectrum in agreement with the observations, and 
with an efficiency of converting the wind kinetic energy in 50-300 keV emission of order 1%. For the same 
set of model features the interval time between peaks and pulse fluences have distributions consistent with the 
log-normal distribution observed in real bursts. 
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1. INTRODUCTION 

The Gamma-Ray Burst (GRB) light-curves are complex and 
irregular, without any systematic temporal features (Fishman & 
Meegan 1995) and an understanding of the origin of the tempo- 
ral behavior of GRBs remains an open issue. Statistical studies 
are necessary in order to identify the physical properties of the 
emission mechanism existent in all or a group of GRBs. Re- 
cently Beloborodov et al. 1998, hereafter BSS98, have used 
the Fourier analysis of a sample of long GRB light-curves to 
study the statistical properties of their power density spectra 
(PDS). The PDS features together with other temporal proper- 
ties of the observed GRBs, such as the distributions of the time 
interval between peaks and of the pulse fluence (McBreen et 
al. 1994, Li & Fenimore 1996), can be used to constrain the 
physical characteristics of the GRB source. 

In the framework of the internal shock model, the rapid vari- 
ability and complexity of the GRB light-curves is due to the 
emission from multiple shocks in a relativistic wind (Rees & 
Meszaros 1994, Kobayashi et al. 1997, Daigne & Mochkovitch 
1998). The ejecta are released by the source during a time com- 
parable to the observed burst duration. The instability of the 
wind leads to shocks which convert a fraction of the bulk kinetic 
energy in internal energy at a distance B ~ 10 12 — 10 14 cmfrom 
the central engine. A turbulent magnetic field is generated and 
electrons are shock-accelerated, leading to synchrotron emis- 
sion and inverse Compton scatterings. Within the framework 
of the internal shock model an alternative hypothesis about the 
particle acceleration and radiation emission is the quasi-thermal 
Comptonization proposed by Ghisellini & Celotti (1999), in 
which particles are re-accelerated for all the duration of the col- 
lision. 

In this paper we analyze the features of the GRB light-curves 
arising from internal shock model, in order to identify the pa- 
rameters that affect most strongly the GRB emission (§2). By 
comparing the features of the simulated bursts with the ob- 
served burst PDS and the distributions of the interval time be- 
tween peaks and of the pulse fluence, we constrain some of the 
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physical properties of the ejecta. 

2. OUTLINE OF THE MODEL 

We simulate GRB light-curves by adding pulses radiated in 
a series of internal shocks that occur in a transient, unstable 
relativistic wind. As we showed in PSM99 the observed burst 
variability time-scale depends mostly on the wind dynamics, 
its optical thickness and its radiative efficiency in the BATSE 
window. Here we model the wind dynamics and the emission 
processes as in PSM99, but we include a more accurate treat- 
ment of the photon down-scattering on the cold electrons in the 
wind. We calculate the effect of the photon diffusion through 
the colliding shells and the wind on the pulse duration and on 
the energy of the emergent photon, rather than just attenuating 
the pulse fluence according to the optical thickness of the wind 
through which it propagates. However the contribution of these 
photons to the duration of the received pulses may be impor- 
tant for bursts that are not very optically thin, and the photon 
down-scattering should be taken into account for more reliable 
calculations of GRB light-curves. 

As described in PSM99, the wind is discretized as a sequence 
of N — t w /t v shells, where t w is the duration time of the wind 
ejection from the central source and t v «t w is the average in- 
terval between consecutive ejections. The shell Lorentz factors 
Ti are random between T m and Tm, where Tm can be con- 
stant during t w ("uniform wind") or modulated on time scale 

t w ("modulated wind"). The shell mass Mj is drawn from 
a log-normal distribution with an average value M — M w /N 
and a dispersion <j m = M, where M w is the total mass in 
the wind, allowing thus the occasional ejection of very mas- 
sive shells. The total mass is determined by requiring that 
Y^iLi MiTiC 2 — L w t w where L w is the wind luminosity. The 
time interval Ai,; between two consecutive ejections i and i + 1 
is proportional to the z-th shell energy, resulting in a wind lu- 
minosity constant throughout the entire wind, and equal to a 
pre-set value L w . This implies that more energetic shells are 
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followed by longer "quiet" times, during which the "central en- 
gine" replenishes. 

Given the wind ejection, we calculate the radii where internal 
collisions take place and determine the emission features for 
each pulse: observer frame duration, fiuence, and photon arrival 
time T b, accounting for relativistic and cosmological effects. 
The peak photon flux for each pulse is calculated assuming the 
pulse shape that Norris et al. (1996) identified in the real bursts, 
described by a two-sided exponential function. The addition of 
all the pulses, as seen by the observer in the 50-300 keV range 
(the 2 nd and 3 rd BATSE channels) gives the burst 7-ray light- 
curve, that is binned on time-scale of 64 ms and is used for the 
computation of the power spectrum. 

For each collision there is a reverse (RS) and a forward 
shock (FS). The shock jump equations allow the calculation 
of the physical parameters of the shocked fluids (Panaitescu & 
Meszaros 1999), determine the velocity of these shocks v s h, the 
compression ratio, the thickness A of the merged shell at the 
end of the collision, and the internal energy in the shocked fluid 
U' s (primed quantities are measured in the co-moving frame). 
The accelerated electrons - a fraction £ e of the total number - 
have a power-law distribution of index — p, starting from a low 
Lorentz factor j m . Assuming that the energy stored in elec- 
trons is a fraction e e of the internal energy, we calculate j m (see 
PSM99). The magnetic field B is parameterized through the 
fraction es of the internal energy it contains: B 2 = 8weBU' s . 
We assume that between two consecutive collisions the thick- 
ness of the shell increases proportionally to the fractional in- 
crease of its radius dA/A oc dR/R. The shell internal en- 
ergy increases in each collision by the fraction of U' s that is not 
radiated, and decreases during the expansion due to adiabatic 
losses. 

The shock-accelerated electrons radiate and the emitted pho- 
tons can be up-scattered on the hot electrons (7 e >> 1) or 
down-scattered by the cold ones (j e w 1). Far from the 
Klein-Nishina regime the optical depth to up-scattering is r ic = 
&ThQen' e miri(ci', A'), where n' e is the co-moving electrons 
density and i! = t' sy /(l+y) is the radiative time scale, with t' sy 
the synchrotron cooling time and y the Comptonization param- 
eter (for r ic < 1, y = j^Tic). The optical thickness r c for the 
cold electrons within the emitting shell is evaluated by taking 
into account the cold electrons within the hot fluid, those that 
were accelerated but have cooled radiatively while the shock 
crossed the shell, and those within the yet un-shocked part of 
the shell. 

A fraction min(l, Tj c ) of the synchrotron photons is inverse 
Compton scattered n ic = max(l,T 2 c ) times, unless the Klein- 
Nishina regime is reached. The energy of the up-scattered pho- 
ton and the ratio of the Compton to synchrotron power can be 
cast in the forms: 
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which take into account the upper limits imposed by the 
Klein-Nishina effect. Figure lc shows the evolution of the 
synchrotron and inverse Compton peak energies during the 
wind expansion: the energy is lower for larger collisions radii, 
due to the increased shell volume and the less relativistic 
shocks, which lead to lower magnetic fields and electron ran- 
dom Lorentz factors. 



The duration STq of the emitted pulse (i.e. ignoring the dif- 
fusion through optically thick shells) is determined by (1) the 
spread in the photons arrival time STg w i?/(2T 2 c) due to 
the geometrical curvature of the emitting shell, (2) the shock 
shell-crossing time 5T& = A/\v s h — i>o|, (where vq is the 
shell pre-shock flow velocity), and (3) the radiative cooling time 
5T 7 w t' /T, which we add in quadrature to determine STq. As 
shown in Figure lb, all these time scales increase on average 
with radius: STg is proportional to R, STa increases due to the 
continuous widening of the shell, and <5T 7 is longer for later col- 
lisions because -f m and B are lower. For ( e = 1, e e ~ 0.25 and 
£b ~ 0.1 the radiative cooling time is negligible respect to STg 
and STa for collisions occurring at R < 5 x 10 14 cm, while 
for larger radii ST^ is the dominant contribution to the pulse 
duration (Figure lb). For the assumed linear shell broadening 
between consecutive collisions we find numerically that the an- 
gular spread and shock-crossing times are comparable during 
the entire wind expansion. 

The optical thickness t c is mainly determined by the wind 
luminosity L w and the range of Lorentz factors in the wind. 
In Figure la for 30 < T < 1000 and L w = lO^ergs" 1 
most collisions occur at R = 5 x 10 13 - 10 15 cm where 
the emitting shells are optically thin. For lower Lorentz fac- 
tors (5 < T < 300) the collisions take place at smaller radii 
(R = 10 12 — 10 13 cm) and the wind is optically thick (Fig- 
ure Id). When r c > 1 photons are down-scattered by the cold 
electrons before they escape the emitting shell, leading to a de- 
crease in the photon energy and an increase of the pulse du- 
ration. For down-scatterings occurring in the Thomson limit 
(e' <C m e c 2 /7e) the energy of the emergent photon can be ap- 
proximated by e' ds = e'(l — e' /m e c 2 ) Tc , where r 2 is the aver- 
age number of scatterings suffered by a photon. For more en- 
ergetic photons, we evaluate e' ds numerically, because the cross 
section depends on the photon energy and changes after each 
photon-electron interaction. For the set of parameters consid- 
ered in this paper, the Thompson limit is usually a good approx- 
imation to treat the down-scattering of the synchrotron photons 
during all the wind expansion. For the smaller collision radii 
the inverse Compton emission peaks at large comoving frame 
energies and the general case has to be considered. Figure 1/ 
shows the evolution of the synchrotron and inverse Compton 
observer frame peak energies for a thick wind. At R w 10 12 
cm, t c 10 3 and the ~ 10 keV synchrotron emission is down- 
scattered by an order of magnitude, while <~ 100 MeV inverse 
Compton radiation is down-scattered to w 10 keV. 

We approximate the increase in the pulse duration due to the 
diffusion through optically thick shells by the time STd it takes 
to a photon to diffuse through them, which we add to STq to de- 
termine the observed pulse duration ST. In the Thompson limit 
STd ~ 5t c A/(24c); in the general case the diffusion time is 
given by 5T d = {J^Zi Wi)]' 1 } A/(2c), where r( Xi ) is the 
optical thickness for the i-th scattering and n s is the number of 
down-scatterings on the cold electrons, evaluated requiring the 
photon random walk equal to the shell width. Figure le shows 
the evolution of the pulse duration during the wind expansion: 
for smaller collision radii STd > STq and the pulse duration 
is determined by STd which decreases with R. For larger radii 
ST d < ST , thus ST w STq and increases with R. 

For a given pulse, we add to the pulse duration the diffusion 
time it takes the photon to propagate through all the shells of 
optical thickness above unity. As shown in Figure Id, the wind 
optical thickness is 1-2 orders of magnitude smaller than the 



3 



optical thickness of the emitting shell (r c ). Nevertheless the 
photon diffusion through the optically thick shells in the wind 
can contribute up to 30% to the pulse duration because of the 
broadening of the shell width during the wind expansion. 

The 30-500 keV pulse energy is a fraction of the kinetic en- 
ergy of the colliding shells, equal to the product of the dynami- 
cal (e^), the radiative (e r ), and the window efficiency (e^). 

1) The dynamical efficiency is the fraction of the kinetic en- 
ergy that is converted to internal, and is given by the energy 
and momentum conservation in the collision of a forward shell 
(Mf,T f ) caught up by a back shell (M b , T b > T f ): 
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is the final Lorentz factor of the merged shell. The Cd decreases 
with T b /T f and is maximized by M b — Mf, so the inner col- 
lisions, for which the difference in the shells Lorentz factor is 
larger, are the most dynamically efficient, with 0.1 . Dur- 
ing the wind expansion the collisions diminish the initial dif- 
ference in the Lorentz factors and the dynamical efficiency de- 
creases to 1 % or less. As show in the next section, a modulation 
in the ejection Lorentz factor is necessary to dynamically effi- 
cient collisions at larger radii. 

2) The radiative efficiency is the fraction of the internal en- 
ergy converted in radiation, and is given by: 
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where t a d ~ R/c is the adiabatic time-scale. The radiative 
efficiency decreases during the wind expansion and it's upper 
limit is the fraction e e of internal energy stored in electrons. 
For magnetic fields not too far from equipartition, the radiative 
timescale is determined by the synchrotron losses. 

3) The window efficiency is the fraction of the radiated energy 
that arrives at observer in the 50-300 keV band. The calculation 
of e w is based on the approximation of the synchrotron spec- 
trum by three power-laws, with breaks at the cooling frequency 
v c and the peak frequency v sy (at which the 7,„-electrons radi- 
ate). If the v c < v sy then the shape of the spectrum is given 
by: 
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where p is the index of the assumed power-law electron distri- 
bution. If v sv < v c then 
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The inverse Compton spectrum has the same shape but is 
shifted to higher energy by the factor implied by equation ([l]). 
For optically thick emitting shells we approximate the burst 
spectrum as given in equations and (fy, using the down- 
scattered cooling and peak frequencies. 



3. EFFECT OF THE PHYSICAL PARAMETERS ON THE GRB POWER 
DENSITY SPECTRUM 

In this section we analyze the effect of the model parameters 
on two distributions that characterize the GRB temporal struc- 
ture: the PDS and the distribution of the interval S p between 
peaks. The relevant model parameters describe the wind ejec- 
tion (t v , t w , T mm , T max and L w ) and the energy release (e e , 
£ e and es). In order to diminish the large PDS fluctuations, 
in this section we use power spectra that are averaged over 10 
peak-normalized bursts. The light-curve peaks are identified 
with the peak finding algorithm (PFA) described by Li & Feni- 
more (1996). For each time bin T p with a photon flux C p higher 
than those of the neighboring time bins, we search for the times 
T\ < T p and T-2 > T p when the photon fluxes C\ and C2 sat- 
isfy C p — Cifl > N vary /Cp~. A peak is identified at T p there 
are no time bins between T\ and T 2 with photon fluxes higher 
than Cp. The light-curve valleys are identified as the minima 
between two consecutive peaks. 

The energy release parameters e e , ( e , and 63 determine the 
50-300 keV radiative efficiency of the pulses. In an optically 
thin wind, the parameter that affects mostly the window effi- 
ciency is the electron injection fraction £ e (for an optically thick 
wind the photons are down-scattered before they escape the 
shells and the window efficiency depends also on r c ). Figures 
2a and 2b show the PDS and the 5 P distributions for a thin wind 
with L = 10 52 ergs" 1 , T m = 30, and T M = 800, and for two 
different Ce (1 and 10 -3 ). In both cases synchrotron emission 
is the dominant radiative process the inverse Compton contri- 
bution to the total emission being 10% for £ e = 1 and less then 
0.1% for ( e = 10~ 3 . For £ e = 1 the synchrotron emission lies 
mainly below the BATSE window (Figure 2c), the window ef- 
ficiency decreases from shorter to longer pulses, the light-curve 
is formed by pulses with a duration of ST w 10~ 2 s, and with 
an average difference in the photon arrival time AT w 0.2 s. 
Because AT ^> ST, the distribution of intervals between peaks 
is determined by the pulses arrival times and peaks at 0. 1 - 0.2 
s. If (e = 10~ 3 the synchrotron emission is above the BATSE 
window for ST < 0.3s and the window efficiency is maximized 
for ST 1=3 0.2 - 0.4 s. The light curve is formed by longer 
pulses, the lower frequency power in the PDS increases and the 
interval time between peaks shifts to longer time-scale. 

The 50-300 keV efficiency of the synchrotron and the in- 
verse Compton emissions is determined by the strength mag- 
netic field B. While for es > 0.1 the emission is dominated by 
synchrotron radiation, for values of the magnetic field well be- 
low equipartition (eg < 0.01) the burst emission is dominated 
by the inverse Compton. Because the shape of the PDS does 
not change, we conclude that the PDS is not sensitive to es and 
the relative contribution of synchrotron and inverse Compton in 
the light-curve. 

The ejection parameters determine the dynamics of the wind 
and the evolution of the pulse dynamical efficiency e<j. The lat- 
ter reflects the evolution of the differences between the Lorentz 
factors of a pair of colliding shells. The first collisions re- 
move the initial random differences, and the merged shells have 
Lorentz factors near the ejection average value T = (T m + 
Tm)/2. If the wind is uniform T is the same for all the shells, 
resulting in a steady decrease of during the wind expansion. 
If the range of shell ejection Lorentz factors is variable on time 
scale of the order of t w (a "modulated" wind), T reflects the 
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initial modulation in Tm and large radii collisions that are dy- 
namically efficient are still possible. 

Figure 3a shows the effect on the PDS of square-sine mod- 
ulations of the upper limit T M with periods P = t w and 
P = t w /A. The i-th shell ejection Lorentz factor is given by 

r, = r m + a 4 sin 2 (j^j (r M -r m ), (8) 

where aj is a random number between and 1, and u> = 
2irt w /P. The modulation shifts the power from high to low 
frequencies, and the magnitude of this shift depends on the 
modulation period. If P = t w the effect of the modulation 
for interaction radii less than rts 10 14 cm (corresponding to 
ST w 1 s) is negligible and the wind evolves as in the uniform 
case: the td decreases from 5% to 0.2% when ST increases 
from 0.01 s to 1 s (Figure 3c). For R <; 10 14 cm the modu- 
lation becomes relevant: the wind is formed of groups of few 
massive shells with different Lorentz factors. The dynamical 
efficiency remains constant for subsequent collisions between 
massive shells, which yield long pulses (ST = 0.3 - 10 s) that 
carry a substantial fraction of the total burst fluence. 

Figure 3d shows that the dependence on ST of the syn- 
chrotron efficiency e sy of the FS pulses has a similar behav- 
ior as that of a, because the internal energy density in the 
shocked plasma depends on e<j. For an higher internal energy, 
the minimum electron Lorentz factor 7 TO increases, leading to 
a higher energy emission and a shorter radiative cooling time- 
scale. Therefore the synchrotron efficiency remains constant 
on the same range of ST where is constant the dynamical ef- 
ficiency, contributing to a shift of power to low frequencies in 
the PDS. 

The optical thickness of the wind depends mostly on the 
range of shell Lorentz factors (T m — Tm) and on the wind lumi- 
nosity (L w ). Figure 4a shows PDSs for two ranges of Lorentz 
factors, 30-1000 and 10-150. In the former case the wind is 
essentially optically thin, and the photon diffusion does not af- 
fect the pulses duration ST, that increases with R (Figure 4c) 
from 0.01 s to 1 s, between 10 13 and 10 15 cm. In the latter case 
the wind is optically thick, in 80% of the collisions r c > 1, 
and the pulse duration is given by the diffusion time: STd de- 
creases from w 5 s to rs 0.6 s between R = 3 x 10 12 cm and 
R = 10 13 cm, where ST is determined mainly by the shell cur- 
vature and thus increases with R. For the optically thick wind 
the long pulses are generated at smaller R where the efficiency 
has the maximum value, and the pulse energy increases with 
ST (Figure Ad). The PDS has more power at low frequency and 
the time intervals between peaks are longer than in the optically 
thin case (Figure 46). The 50-300 keV efficiency is of the same 
order for the two cases: 4 x 10~ 3 and 5 x 10~ 3 for a range of 
Lorentz factors of 30-1000 and 10-150, respectively. 

An increase in the wind luminosity has a similar effect on the 
PDS shape as a decrease in T m and Tm- In the latter case the 
wind becomes thicker because the shells are more massive. 

The variability time scale t v affects the dynamical evolution 
of the shells in the following way. If the time intervals between 
successive ejections delays decreases then the collisions occur 
at smaller radii, where the wind is more optically thick. The dif- 
ferences between the Lorentz factors diminish faster (there are 
more shells for smaller t v ), reducing the dynamical efficiency 
for short pulses. For the modulated wind this effect is more 
relevant than in the random case. The duration t w of the wind 
ejection determines mainly the number of shells, and changes 



in t w do not affect much the evolution of uniform winds. How- 
ever, for a modulated wind, t w also determines the number of 
periods in the Lorentz factor (if the duration of a period is inde- 
pendent of t w ), influencing thus the clumping of shells. 

The burst redshift determines the co-moving energy range 
which is redshifted into the observing range, leading to a 
change in the total pulse efficiency, and altering the observed 
pulse duration. Obviously, by increasing the burst redshift, 
power is shifted from higher to lower frequencies. 



4. COMPARISON WITH THE OBSERVATIONS 

An analysis of the PDS of real bursts was presented by 
BSS98. They calculated the Fourier transform of 214 long 
(Tg > 20 s) and bright burst, and have found that the average 
PDS is a power-law (Pj oc /~ 5 / 3 , / is frequency) over almost 
two orders of magnitude in frequency, between 0.02 Hz and 2 
Hz, where a break is observed, indicating a paucity of pulses 
with duration less than m 0.5 s. The distribution of intervals 
between peaks has been studied by McBreen et al. (1994) and 
by Li & Fenimore (1996), who showed that the distributions of 
the pulse fluence S p and of the time interval S p between peaks 
are consistent with a log-normal distribution. 

As was shown in the previous section, if the wind is optically 
thin and the ejection features are random, the pulse duration 
increases with the collision radius and the emission efficiency 
decreases during the wind expansion. The short inner colli- 
sions yield most of the 50-300 keV burst emission and the in- 
ternal shock model predicts a flat PDS with equal power at low 
and high frequency. Thus, in order to explain the observed be- 
havior, we need a configuration of the parameters which shifts 
power from the short to the long time-scales in the light-curves. 
Moreover, the S p distribution is not log-normal: Figures 26, 3b, 
and 46 show that in GRBs arising from optically thin, uniform 
winds there are too many short intervals between peaks respect 
to a Gaussian log S p distribution. 

PSM99 have identified three possible ways to explain the 
deficit of pulses with ST < 1 s: 

(1) a reduction in the electron injection fraction. This in- 
creases the photon energy, reducing the window efficiency of 
the short pulses (causing the high energy break) and increasing 
that of the longer ones. However the behavior of the PDS at 
lower frequency remains flat (see Figure 2a). 

(2) a modulation of the shell ejection Lorentz factor. This 
allows different configurations for the collisions series and a 
higher dynamical efficiencies for longer pulses (see Figure 3a). 

(3) an increase of the optical thickness of the wind. In this 
case the down-scattering suffered by the photons as they propa- 
gate through the wind increases the pulse duration for the small 
radii collisions, which yield the shorter duration pulses (see 
Figure 4a). 

In Figure 5a we show a simulated light-curve for a square- 
sine modulated wind (with P = t w ) The burst 50-300 keV 
efficiency is 1%, and the 90% of the RS and 80% of the FS 
propagate in optically thick shells. If N var = 0.1 (the free pa- 
rameter of the PFA), we find 22 pulses in the light-curve shown 
in Figure 5a. In order to have more peaks we simulate four 
light-curves with the same injection features and wind parame- 
ters and we calculate the interval between peaks S p (Figure 56) 
and peak fluence S p (Figure 5c) distributions. The distributions 
are similar to a log-normal one, and the choice of N var does 
not affect strongly their shape. 
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In order to compare the PDS of the simulated bursts with 
the observed one, we consider an ensemble of cosmologi- 
cal GRBs. Some authors (Totani, 1997, Wijers, et al. 1998, 
Krumholz et al. 1998, Hogg & Fruchter 1999, Mao & Mo 
1999 ) have used a GRB co-moving rate density proportional 
to the star formation rate. Others (Reichart & Meszaros 1997) 
have employed a power-law GRB density evolution with red- 
shift, which was found by (Bagot et al. 1998) to be consistent 
for z < 2 with their results from population-synthesis com- 
putations of binary neutron stars merger rates. Finally, other 
researchers (Krumholz et al. 1998, Hogg & Fruchter 1999), 
have considered a constant GRB rate density. In this work, 
we use the power-law with redshift GRB density evolution 
n c (z) oc (1 + z) D , mainly as a convenient parameterization. 
An n c (z) proportional to the star formation rate would lead to 
different sets of model parameters (see below), but the differ- 
ences are minor, because the two functions differ substantially 
in shape only for z > 1, where there is a strong decrease of 
the co-moving volume per unit redshift and a smaller chance of 
obtaining a burst that has a 50-300 keV peak photon flux below 
1 7 cm~ 2 s _1 (bursts dimmer than this limit are not included in 
the calculation of the average PDS and intensity distribution). 

Given the rate density evolution, the GRB redshift is chosen 
from a probability distribution 

dV n c (z) dV 

t~ x ~r, — ~t~ ' ( 9 ) 

dz 1 + z dz 

where dV/ dz is the cosmological co-moving volume per unit 
redshift 

dV = f c\ 3 [ g0 z-(l- g0 )(V2g z + l-l)] 2 
dz n \H ) q^qoz + iy/^l + zf 

(10) 

We assume qo = 0.5 and Ho = 75 km s 1 Mpc 1 . 

The inferred isotropic 50-300 keV luminosities of the GRBs 
that have measured redshifts span more than one order of mag- 
nitude, therefore the standard candle approximation is not a 
good approximation. We use an un-evolving power-law dis- 
tribution for the wind luminosity: 

oc L^ 13 , L m < L < Lu i (11) 

and zero otherwise. Note that this not the same as assum- 
ing that GRBs have a power-law distribution of their 50- 
300 keV luminosities, as it is usually done (e.g. Reichart & 
Meszaros (1997), Krumholz et al. 1998, Mao & Mo 1999), as 
the relationship between the wind and the 50-300 keV lumi- 
nosities is set by the window efficiency (at the source) and, in 
the case of winds that are optically thick, by the wind optical 
thickness, both of which are dependent on the wind luminosity. 

In finding model parameters that yield bursts consistent with 
the observations, we held constant t v — 25 ms, Lm/Liu = 
100, (3 = 2, e e = 0.25, e B = 0.1, andp = 2.5 . The cho- 
sen t v is short enough to ensure that the observed 2 Hz PDS 
break frequency is below [(1 + z)^] -1 (bursts with z > 3 are 
rarely brighter than 1 7 cm _2 s _1 ), corresponding to the pulses 
that are partly suppressed by the choice of t v . the PDS fre- 
quencies affected by the assumed t v = 25 ms are <; 10 Hz. 
This may suggest a possible explanation for the PDS break 
observed by BSS98: the lack of pulses shorter than ~ 1 s is 
due to the existence of a minimum wind variability time-scale 
of the same order. However, such i„'s would be much larger 



than the dynamical time-scales of plausible GRB progenitors 
(Meszaros et al 1999), and we do not consider this a viable 
possibility. The choices of Lm / L m and (3 are consistent with 
the values found by Reichart & Meszaros (1997), Mao & Mo 
(1999), and Krumholz et al. (1998) from fits to the observed 
intensity distribution. The values chosen for e e and e b are not 
too far from those determined by Wijers & Galama (1999) from 
the emission features of the afterglows of GRB 970508 and 
971214. The above value of electron index p is close to the 
values implied by the observed slopes of the afterglow optical 
decays. 

Figure 6b shows a burst-averaged PDS whose features are 
similar to that found by BSS98 in real bursts. The wind ejec- 
tion is modulated by a square sine (eq. [^]) with a random pe- 
riod between t w /4 and t w . About 40% of the 300 simulated 
bursts have peak photon fluxes brighter than l7cm _2 s . 
Taking into account that the average redshift for these bursts 
is ~z = 0.90 the average burst duration Tt « 1.5(1 + z)t w 
is close to the value T\, = 80 s of the bursts used by BSS98 
(the factor 1 .5 was determined numerically and represents the 
ratio between the burst duration at the source redshift and t w ). 
As can be seen in Figure 6b, P f oc f- 5 /' 3 between 0.04 Hz 
and 2 Hz and falls off steeper at frequencies larger than 2 Hz. 
The model parameters that led to the PDS of Figure 6b yield 
bursts whose integral intensity distribution is shown in Figure 
6a, which consistent with the distribution found by Pendleton 
et al. (1996): excluding the bursts dimmer than 1 7cm~ 2 s -1 , 
the model has \ 2 — 9.5 for 9 degrees of freedom. 

5. CONCLUSION 

We have calculated power density spectra of GRBs arising 
from internal shocks in an unsteady relativistic wind. By study- 
ing how the features of these spectra depend on the model pa- 
rameters (Figures 2, 3, and 4), we have identified a set parame- 
ters (Figure 6) that leads to bursts whose average PDS exhibits 
an f~ 5 ' 3 behavior (where / is frequency) for 0.04 Hz < / < 
2 Hz, as found by BSS98 in real GRBs. Moreover, the integral 
intensity distribution of the simulated bursts is consistent with 
that observed by Pendleton et al. (1996), and the distributions 
of the time intervals between peaks and of the pulse fluences 
are consistent with the log-normal distributions identified by Li 
& Fenimore (1996) in real bursts. 

The characteristics of the modeled bursts with the above 
mentioned features are: (1) a sub-unity electron injection frac- 
tion, required to increase the radiative efficiency of the larger 
collision radii, (2) a modulated Lorentz factor of the ejected 
shells, necessary to increase the dynamical wind efficiency dur- 
ing the wind expansion and, (3) a shells optical thickness to 
scattering on cold electrons above unity, required to increase 
the duration of the pulses as they propagate through the collid- 
ing shells and the wind. 

In the internal shock model, the most efficient collisions, 
with a dynamical efficiency of 10-20% and a radiative effi- 
ciency of 10 — 30%, happen in the first part of the wind ex- 
pansion where the wind optically thickness is higher and the 
angular spread time, the shell shock-crossing time and the elec- 
trons cooling time are shorter (<C 0.5 s). In order to reproduce 
the observed break at 2 Hz in the PDS, we have previously (see 
PSM99) attenuated the fluence of these short pulses according 
to an high wind optically thickness, with a resulting low burst 
efficiency (10~ 4 for an uniform wind and 10~ 3 for a modu- 
lated one). The study of the photon diffusion, presented here, 
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allowed us to find model parameters that yield an 1 % efficiency 
of converting the wind kinetic energy into 50-300 keV emis- 
sion. For an optically thick wind, the pulse duration of the first, 
efficient collisions at small radii is determinated by the time 
the photons take to escape the shells, that depends only on the 
colliding shells width and optically thickness r c . If r c 3> 1 
the diffusion time for the efficient collisions is > 0.5 s and the 



simulated average PDS shows the break at 2 Hz with a burst 
efficiency close to the maximal value (few %) admitted by the 
model (see also Kumar 1999). 

This research is supported by NASA NAG5-2857, NSF 
PAY94-07194 and the CNR. We are grateful to Martin Rees, 
Stein Sigurdsson and Marco Salvati for stimulating comments. 
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Fig. 1. — Upper panels: optically thin wind with L w = 10 53 ergs -1 and 30 < T < 1000. Panel a shows the dependence on the collision radius R of the 
optical thickness t c for scattering on the cold electrons inside the emitting shell, and that of the rest of the wind (t w ). Panel b illustrates the i?-dependence of the 
pulse duration and of the terms that contribute to it; while graph c shows the R-dependence of the synchrotron and inverse Compton peak energies. Lower panels: 
optically thick wind L w = 10 53 ergs -1 and 5 < V < 300. Graph d shows the R-dependence of r c and t w , graph e shows the evolution of the pulse duration 
ST = &Td + (<5Tg + + (5T?) 1 / 2 and the contribution of the diffusion time through the wind <5Tj (i.e. excluding the emitting shell). Graph / shows the 
down-scattered synchrotron and inverse Compton energy peaks versus the collision radius. The dashed lines in panel / show the evolution of the synchrotron and 
inverse Compton peak before the down- scattering, and the dotted lines in panel c and / show the BATSE window. Parameters: t v = 0.02 s, t w = 20 s, e e = 0.25, 
eg = 0.1, fe = 1, burst redshift z = 1. Only a small fraction of the total number of pulses is shown; the density of the points illustrates the radius distribution. The 
curves shown are log-log space fits for the most efficient pulses. The actual values are scattered around the fit. 
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Fig. 2. — Effect of the injection fraction f e on the burst PDSs (panel a), distribution of interval between peaks (panel b) for Nt ~ 130 peaks, energy of the 
synchrotron spectrum peak (panel c), and window efficiency (panel d). Parameters: L w = 10 52 erg s _1 , T m = 30, T M = 800, other parameters are as for Figure 
1. Panels c and d show log-log space polynomial fits, illustrating thus only the trends. The RS emission is represented with circles and triangles, while the FS one is 
shown with crosses and stars. 
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Fig. 3. — Effect of a modulation of the wind ejection. Panels a and b show a comparison between the PDSs and the S p distributions for an uniform wind (solid 
line), and a square-sine modulated ejection with period t w (dotted line) and t w /4 (dashed line). Parameters: L w = 10 erg s , 1 = 30, T M = 1000, other 
parameters are as for Figure 2. Lower panels: the dependence of the dynamical efficiency (panel c), and FS synchrotron efficiency (panel d) versus the pulse duration. 
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Fig. 4. — Effect of the optical thickness. Panels a and b show a comparison between the PDSs and the distributions of S p for an optically thin wind with 
30 < T < 1000 (Nt = 70 peaks), and an optically thick one with 10 < T < 150 (Nt = 40 peaks). L w = 10 53 ergs -1 , other parameters are as for Figures 
2 and 3. The dependence of ST on the collision radius and of the pulse 50-300 keV fluence on ST are shown in panels c and d respectively, for the optically thick 
(triangle) and thin (circle) winds. The pulse fluence is the fraction of the total fluence in the light curve carried by each pulse. 
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Fig. 5. — Top panel: light-curve for a modulated wind with P = t w and T m = 5, Tm = 150, L w = 2 X 10 52 ergs" 1 , t v = 25 ms, t w = 30 s, z = 1, 
£e = 0.1, e e = 0.25, eg = 0.1 . The photon flux is normalized to its highest value. The triangle the peaks (38) and the valleys selected with the PFA (N var = 0.1). 
Bottom: the distributions of the interval 8 P between peaks (b) and of the peak fluence S p (c), the latter representing the fraction of the total fluence between two 
consecutive valleys. The distributions are calculated for four bursts with identical ejection features and wind parameters, total numbers of peaks is 80. The log-normal 
distributions that fit best the numerical results are shown with dotted curves. 
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U<)N/(d<)N 

Fig. 6. — Left panel: intensity distribution of simulated bursts (open squares) compared to that of the observed ones, taken from Pendleton et al. (1996) (filled 
circles). Both distributions are normalized to the number of bursts brighter than 1 7 cm _2 s _1 . Right panel: the averaged PDS of the simulated bursts, compared 
with the shape of the PDS (thick line) determined by BSS98 for real bursts. The average is done over 112 bursts with a square-sine modulated wind of random period 
between t w /4 and t w . Parameters: L m = 2 X 10 52 ergs" 1 , L M = W0L m , P = 2, D = 3, t w = 30 s, t„ = 25 ms, e e = 0.25, e B = 0.1, £ e = 0.1, Y m = 5, 
T M = 150, D = 3. 



